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ABSTRACT 

We analyze the consequences of models of structure formation for higher-order (n-point) galaxy 
correlation functions in the mildly non-linear regime. Several variations of the standard fi = 1 cold 
dark matter model with scale-invariant primordial perturbations have recently been introduced to obtain 
more power on large scales, R p ~ 20 h~ l Mpc, e.g., low-matter-density (non-zero cosmological constant) 
models, ‘tilted 1 primordial spectra, and scenarios with a mixture of cold and hot dark matter. They 
also include models with an effective scale-dependent bias, such as the cooperative galaxy formation 
scenario of Bower, etal. (1993). We show that higher-order (n- point) galaxy correlation functions can 
provide a useful test of such models and can discriminate between models with true large-scale power in 
the density field and those where the galaxy power arises from scale-dependent bias: a bias with rapid 
scale-dependence leads to a dramatic decrease of the hierarchical amplitudes Qj at large scales, r > R p . 
Current observational constraints on the three-point amplitudes Q$ and S 3 can place limits on the bias 
parameter(s) and appear to disfavor, but not yet rule out, the hypothesis that scale-dependent bias is 
responsible for the extra power observed on large scales. 

Subject Headings: Large-scale structure of the universe — galaxies: clustering 
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1 Introduction 


Recent observations of galaxy clustering in both photometric and spectroscopic surveys have found more 
relative power on large scales, R p ~ 20 h~ l Mpc (h = H 0 / 100 km/sec/Mpc), than that expected in the 
standard cold dark matter (CDM) model of structure formation (e.g., Maddox, etal. 1990, Efstathiou, 
etal. 1990, Baumgart and Fry 1991, Gramann and Einasto 1991, Hamilton, etal. 1991, Peacock and 
Nicholson 1991, Saunders, etal. 1991, Loveday, etal. 1992, Fisher, etal. 1992, Park, etal. 1992, Vogeley, 
etal. 1992, Feldman, etal. 1993). More precisely, the shape of the observed galaxy power spectrum 
P g (k) or of its Fourier transform, the two-point galaxy correlation function £ 9 (r), differs on these scales 
from the standard CDM model prediction. 

Recall that in the standard CDM model, the Universe is spatially flat, with a density Q c dm = 

1 — Qb — 0 95 in non-baryonic, weakly interacting particles which have negligible free-streaming length, 
and the Hubble parameter h = 0.5. Additionally, one posits that the density perturbations responsible 
for large-scale structure are adiabatic and Gaussian, with a scale-invariant primordial power spectrum 
P(k) = (\S k (ti)\ 2 ) ~ k, as expected in canonical inflation scenarios. The present spectrum is related 
to the primordial one through the transfer function, T{k\Ui,h), which encodes the scale-dependence of 
the linear growth of perturbations, (|M*o)| 2 ) = 7 ’ 2 ( fc )(l^( t *)l 2 )- Finally, the galaxy power spectrum is 
related to the density spectrum by a bias factor b g , 

P g (k) = b 2 g T 2 (k)Mti)\ 2 . (1) 

A number of alternatives have been suggested to remedy the shape of the CDM galaxy spectrum, 
each of which involve modifications of one or more of the standard ingredients of the CDM model in 
equation (1). These include models with a lower density of cold dark matter, fi c dmh — 0.2, plus a 
cosmological constant to retain spatial flatness (Efstathiou, Sutherland, and Maddox 1990), and models 
with a mixture of cold and hot dark matter, fl c dm — 0.7, £2hdm — 0.3 (e.g., Schaefer, etal. 1989, 
van Dalen and Schaefer 1992, Taylor and Rowan- Robinson 1992, Davis, etal. 1992, Pogosyan and 
Starobinsky 1992, Klypin, etal. 1992). In these two cases, the transfer function T(k) is flattened on 
scales it -1 ~ Rp compared to standard CDM. CDM models with ‘tilted’ non-scale-invariant, power-law 
primordial spectra, (|£*(i,)| 2 ) ~ k n with n < 1, which arise naturally in several models of inflation, 
have also been recently explored (Adams, etal. 1993, Cen, etal. 1992, Gelb, etal. 1993, Liddle and 
Lyth 1992, Liddle, etal. 1992, Vittorio, etal. 1988). In addition, there is a growing literature on 
models with non-Gaussian initial fluctuations; in some cases, initial skewness and/or kurtosis can lead 
to enhanced structure on large scales (e.g., Moscardini, etal. 1993 and references therein). While such 
models can display interesting behavior of the higher order moments, in this paper we will focus on 
initially Gaussiam fluctuations. 

In all these variations on the CDM theme, one important assumption is left unchanged: that the 
observable galaxy distribution is related through a simple bias mechanism to the underlying matter 
distribution predicted by theory (e.g., Bardeen, etal. 1986). In essence, following Kaiser (1984a, b) 
and Baurdeen (1984), one assumes that galajdes form from peaiks above some global threshold in the 
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smoothed linear density field. In the limit of high threshold and small variance, this model is well 
approximated by the commonly employed linear bias scheme, in which the galaxy and mass density 
fields, 6 g (x) = (n g (x) - n g )(n g and S(x) = (p(x) — p)fp, are linearly related through a constant bias 
factor, 

M*) = M( x ) ■ (2) 

This relation, implicitly assumed in equation ( 1 ), embodies the standard model for biased galaxy for- 
mation. 

Early numerical evidence for biasing came from the CDM simulations of White, etal. (1987), which 
showed that dark matter halos are more strongly clustered than, and thus ‘naturally' biased with re- 
spect to, the mass. However, since galaxy formation is a complex, non-linear process involving both 
gravitational and non-gravitation&l interactions, the relation between the mass and the galaxy distribu- 
tions may be more complicated than in the peak bias model. Even purely gravitational high-resolution 
N-body simulations suggest that virialized halos are not always well identified with peaks in the linear 
density field (Katz, Quinn, and Gelb 1992). 

It is therefore of interest to ask whether a more or less well-motivated modification of the standard 
bias scheme can generate the excess large-scale power within the context of the standard CDM model. 
This idea has been recently studied by Babul and White (1991) and by Bower, etal. (1993) (for 
precursors, see Rees 1985, Silk 1985 and Dekel and Rees 1987). The common thread in these ideas is 
that the bias mechanism can be modulated by environment-dependent effects. For example, in their 
cooperative galaxy formation scenario, Bower, etal. (1993) (hereafter BCFW) suggest that the threshold 
above which perturbations actually form bright galaxies may be lower in large-scale, high-density regions 
than elsewhere. Or perhaps baryons may be inhibited from cooling in regions photoionized by an early 
generation of quasars (Babul and White 1991). The net result of these feedback mechanisms is that the 
transformation from the density field 6(x) to the galaxy field S g (x) becomes non-local (by contrast with 
equation (2)), and the effective bias factor becomes scale-dependent. If the bias factor increases with 
scale, the galaxy spectrum will have more power at large scales, as desired. This modification of the 
standard CDM scenario is fundamentally different from those mentioned above: with scale-dependent 
bias, the extra large-scale power relative to standard CDM is only apparent, in the sense that it is only 
a property of the galaxy field, not the underlying mass density field; by contrast, in the other CDM 
variants (non-zero A, tilt, or mixed dark matter), there is genuine extra power in the density field. 

In this paper, we consider how the higher order irreducible moments of the galaxy distribution can 
be used as a test of models for large-scale structure. We consider the standard CDM model and its 
variants with extra large-scale power (in particular, Qh — 0.2 CDM), as well as a generalized version 
of the non-local, scale-dependent bias scheme embodied in the cooperative galaxy formation (hereafter, 
CGF) model of BCFW in the context of otherwise-standard CDM. Using the results of second-order 
perturbation theory (Fry 1984), we compare in detail the predictions of these models for the three- 
point function £3 with data from the Center for Astrophysics (CfA, Huchra, etal. 1983), Southern Sky 
(SSRS, Da Costa, etal. 1991), and Perseus- Pisces (Haynes and Giovanelli 1988) redshift surveys in the 
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mildly non-linear regime (£2 < !)■ Since £3 is of second-order in the density perturbation amplitude 
for initially Gaussian fluctuations, for self-consistency we must generalize the models to include the 
possibility of non-linear (as well as non-local) bias and extend them from Gaussian to hierarchical 
matter fields. [We will use the well known result that, at least in the mildly non-linear regime, the matter 
field evolved gravitationally from Gaussian initial conditions leads to hierarchical statistics of the form 
(6 J ) oc ((5 2 )* 7 ” 1 (cf. Fry 1984, Goroff et al 1986, Bemardeau 1992).] The allowance for non-linear bias 
introduces an additional dimensionless parameter into the model. Even with this additional degree of 
freedom, we find that the CGF model tends to require rather large values of the bias parameter in order 
to match the 3-point function data, because scale-dependent bias modifies the correlation hierarchy, 
leading to a dramatic decrease of the hierarchical amplitudes Qj at large scales, r > Rp . In the context 
of standard CDM, such a high bias is in conflict with the COBE DMR observations of microwave 
anisotropy on large scales. We show that observations of the 3-point function in Fourier space, Q(k), on 
the largest scales accessible to current redshift surveys should provide a definitive test of the CGF model 
and of more general models with scale-dependent bias. Our basic conclusion is that the scale-dependent 
bias solution to the problem of extra large-scale power affects the 3-point functions very differently from 
models with genuine extra power (such as CDM with Qh = 0.2). Thus, the higher-order correlations 
provide an important test to distinguish between different solutions of the extra power problem. 

The paper is organized as follows. In section II, since it may be less familiar to the reader, we 
briefly review and generalize the CGF model and recapitulate the results of BCFW, demonstrating 
the enhancement of the two-point function on large scales required to fit the APM angular correlation 
function data. In section III, we review the results on the 3-point and higher order correlations in 
perturbation theory, focusing on the evolution of an initial Gaussian density field into a hierarchical 
field. In section IV, we study the higher order moments in the CGF model. Self-consistency demands 
that we further extend the model to include non-linear bias. In section V, we compare the standard 
CDM, low-density CDM, and CGF-modified CDM predictions to the data on the 3-point function from 
the CfA, SSRS, and Perseus- Pisces redshift surveys and we conclude in section VI. 

2 Cooperative Galaxy Formation and Scale-Dependent Bias 

The cooperative galaxy formation (CGF) model of BCFW is a simple phenomenological prescription 
for obtaining a scale-dependent bias. It starts with the standard assumptions of the CDM model, but 
the biasing mechanism is modified from the high peak threshold scenario. In the standard peak bias 
model (Kaiser 1984a, Bardeen, etal. 1986), the sites of galaxy formation are identified with peaks 
of the smoothed linear density field. That is, one convolves the initial density field with a filter of 
characteristic scale R g ~ 1 h~ l Mpc, and then identifies galaxies with peaks of the smoothed field above 
some threshold vo y i.e., with density maxima satisfying 6(® p *) > vo, where o\* = (( p — p) 2 )/p 2 is the 
variance of the smoothed field, and v sets the threshold height. (Hereafter, we implicitly assume the 
field 6(x) is smoothed on the scale R g .) For example, for an infinitely sharp threshold, the galaxy field 
is 6 g (x p k) = 0(6(x p k) - 1 / 0 ). The combination of the threshold peak height v and the spatial smoothing 
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scale R g is chosen so that the density of peaks reproduces the observed abundance of luminous galaxies; 
moreover, these parameters are taken to be global, spatially invariant quantities. In the limit of high 
threshold (v » 1) and small variance, the two- point correlation function of the peaks is enhanced over 
that of the mass by an approximately constant factor (Kaiser 1984a), 


Zpk(r;v) 



( 3 ) 


where £(r) = {6(x)6(x + r)). Since £(r) is quadratic in the density field, this is equivalent to the linear 
bias model of equation (2), with the identification of the bias factor as b g = {vjo). Following Kaiser 
(1984a) and BCFW, we will apply this model to regions above the threshold, S(x) > va , rather than 
to maxima; this simplifies the model while retaining its important features. 

BCFW extend the standard bias model by replacing the universal threshold v with a threshold that 
depends on the mean mass density in a surrounding ‘domain of influence’ of characteristic size R s > R g . 
The motivation is to model the possibility that peaks form galaxies more easily (or perhaps form brighter 
galaxies which are included in a magnitude-limited catalog) if there are other peaks nearby-thus the 
name cooperative galaxy formation. Specifically, they assume that galaxies form from regions satisfying 


6{x) > va - kS(x;R s ) , ( 4 ) 

where £(ar; R s ) is the density field smoothed on the scale R s , and k is the modulation coefficient of the 
threshold. If k > 0, the threshold for galaxy formation is lower in “protosupercluster” regions than in 
“protovoids” . The parameters R s and k parametrize the scale and strength of cooperative effects; they 
are also constrained by the observed galaxy abundance. 

The model of equation (4) is equivalent to applying the standard threshold bias model to the new 
density field defined by 

6'(x) = 6(x) + «£(*; R a ) , (5) 

that is, to imposing the condition 8 f > vo , Note that 8 f is a Gaussian random field if the underlying 
density field 8 is Gaussian. Here we consider a generalization of the CGF model: instead of applying a 
sharp threshold clipping to 8*(x), we assume that the galaxy field is an arbitrary continuous function 
of the field 6 * , 

M*) = /(«'(*)) = / [«(*) + *6(x; *,)] . (6) 

For example, in the limit of high threshold, for the standard bias model the function / is approximately 
an exponential, f(x) = exp(vx/o) (Kaiser 1984b, Politzer and Wise 1984). We assume that / is 
expandable in a Taylor series in its argument, 

6, = m - E I 4 '* • (?) 

fc=i *• 

BCFW compute the two-point correlation function for the CGF model on large scales, using the CDM 
density spectrum derived from linear perturbation theory. In this regime, our generalized CGF model 
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reduces to the linear bias model applied to the field 6 \ 

6g(x) = b g 6'(x) = bg [<5(x) + k6(x; i? s )j (8) 

where we have identified b g = b\. That is, by working only to linear order in perturbation theory, one 
should self-consistently include only the first (linear) term in the series of equation (7). Conversely, when 
we consider second order perturbations below, we can and should include the possibility of quadratic 
(k = 2) bias. 

Comparing equation (8) with equation (2), it is clear that cooperative effects boost the galaxy power 
spectrum on large scales relative to the standard global bias model. Taking a Gaussian filter for the 
smoothed density field, 

6 (*; R a ) = (2t *R *) " 3/2 J d 3 rS(r)ex p . ( 9 ) 

the Fourier transforms of the density fields satisfy 

6'(k) = S(k) [1 + /c£(fc)] , (10) 

where Q is the Fourier transform of the window filter in 6(x; R a ), 

Q{k) = G(k) = e~( kR ‘ )2 1 2 , (11) 

with k = \k\. The galaxy power spectrum, P g (k) = (|6 fl (fc)| 2 ), is thus related to the density power 
spectrum, P(k) = (|6(fc)| 2 ), by 

Pg(k)=b 2 g P , (k) = b 2 g[l + KG(k)] 2 P(k) = bl fr (k) P(k) . (12) 

This expression makes manifest how cooperative effects result in an effective scale-dependent bias, 
i> e ff(fc) = b g [ 1 + KQ(k)}. On small lengthscales, /c" 1 < R s , equation (12) implies the usual bias factor, 
&eff(k — ► oo) csf b g , while on large scales, A: -1 > R a , the effective bias factor is increased to 6 e ff(A: — ► 0) ~ 
b g (l + ^c). In the parameter range studied by BCFW, the choice k = 2.29, R $ = 20/i _1 Mpc appears to 
give the best fit to the observed extra large-scale power for CDM when compared to the APM angular 
correlation function, and we shall focus mainly on this case. We see that this choice boosts the galaxy 
power spectrum on scales k < 0.05/i Mpc -1 by over a factor of ten. 

To see what these effects look like graphically for the CDM model, we consider the linear CDM 
density power spectrum of Davis, etal. (1985), 

P(k) = Ao\k (l + + Jilp + jjip) , (13) 

where the wavenumber k is in units of h Mpc -1 . Here the normalization is set as usual in terms of the 

1 /2 

variance of the linear mass fluctuation within spheres of radius 8 h -1 Mpc, erg = ((6M/M g^-iMpc’ 
where 

oo 

<r% = ^2 Jdkk 2 P(k)W 2 (kR) , (14) 

o 
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and the top-hat window function 


W(kR) = 


3 

(kR) 3 


filters out the contribution from small scales. 


(sin kR — kR cos kR) (15) 

For standard CDM with 0 h = 0.5, this gives A = 


2.76 x 10 5 ( h~ l Mpc) 3 . 

Substituting the CDM power spectrum with Q/i = 0.5 into equation (12), we find the galaxy two- 
point correlation function for the CGF model 


. , . 1 f . o sin kr _ . 

^( r ) “ 2^2 J dkk ^r p s( k ) ’ (i®) 

shown in Fig. 1 (the curve labelled CGF, with k = 2.29, R 9 = 20 h~ l Mpc). Note that we actually 
plot ^(r)/(6^<7g) 2 , where b g is the constant factor in equation (8). Redshift surveys of optically selected 
galaxies (in particular the CfA and Stromlo-APM surveys) indicate that the variance in galaxy counts 
on 8 /i -1 Mpc scale is of order unity. Thus, in a linear, scale-independent bias model, the bias factor for 
these galaxies would be expected to be — 1 / erg; for other galaxy populations, however, b ga ic rg may 
differ from unity. For comparison, in Fig. 1 we also show the two-point function for standard CDM 
(Qh = 0.5, k = 0) and for a low-matter-density CDM model (fift = 0.2). Both the CGF model and 
the low-density CDM model have sufficient relative large-scale power to approximately reproduce the 
observed galaxy angular correlation function w(0) inferred from the APM survey (BCFW, Maddox, etal. 
1990, Efstathiou, Sutherland, and Maddox 1990). This level of extra power is also broadly consistent 
with that inferred from the power spectrum of IRAS galaxies (Feldman, etal. 1993, Fisher, etal. 1992) 
and the redshift-space two-point function £(s) inferred from the Stromlo-APM survey (Loveday, etal. 
1992). The CGF curve in Fig. 1 should be compared to that in Fig. 2 of BCFW. Note that the linear 
bias approximation used here (equation (8)) differs from the non-linear threshold formula of BCFW 
(Cf. their eqn.(10)), but that our final result for £(r) is very similar to theirs. 

Thus, cooperative effects can mimic extra large-scale power in the galaxy two-point function, while 
the other remedies for CDM, such as low-density, mixed dark matter, or tilted (n < 1) models, have 
genuine extra large-scale power in the spectrum. How can we discriminate between these choices for 
extra large-scale power, that is, between real power and the illusion of power? Below, we argue that 
the three-point function can provide a distinguishing test, at least for models with Gaussian initial 
fluctuations. The reason is that the galaxy three-point function induced by gravitational evolution 
depends in large measure on the two-point function of the mass . 

Before turning to higher order correlations, we remark that the treatment given here and below 
applies more generally than to the CGF model, and in fact to any model with scale-dependent bias. 
That is, the chain of reasoning above is invertible: if the galaxy and density power spectra are related 
by a scale-dependent bias, P g (k) = b 2 (k)P(k), we can always think of the galaxy field 6 g (z) as arising 
from some non-local transformation of the density field 6(x). To see this, let b 2 (k) = b 2 f 2 {k), where 
b g is a constant and we assume that f 2 {k) has a limit, f 2 (k — ► oo) = 1. Then consider the field 
S f (k) = 6 g (k)/b g = f{k)6(k ). We can write f(k) = 1 + *G(fc), where lim G(k — > oo) = 0, and we can 
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choose k such that lim G{k -> 0) = 1, so that 6'(k) = 6 (fc)[l + kG(A:)]. Comparing with equation (10), 
we see that this expression, where G(k) is interpreted as the Fourier transform of some window function 
(which in general will not be a Gaussian), is all that we need for the results discussed here and below 
to go through. Provided the function b(k ) is not too pathological, this Fourier transform should exist. 


3 3-point correlation function in perturbation theory 

We want to consider how scale-dependent bias, as embodied for example in the CGF model, affects 
the higher order correlation functions. The motivation for this study is that the galaxy three-point 
function is observed to scale in a particular way with the two-point function, and both perturbation 
theory and N-body simulations show that this scaling can arise via non-linear gravitational evolution 
from Gaussian initial fluctuations. Since scale-dependent bias introduces a different scale behavior into 
the problem, we would expect it to be manifest as a change in the scaling behavior of the higher order 
correlations. We will work in the context of second-order perturbation theory (Fry 1984), the results of 
which we review here before discussing how they are modified by scale-dependent bias. The perturbative 
approach should be valid in the mildly non-linear regime, S < 1. In the range where they overlap, the 
second-order perturbation theory results below for 53 in standard CDM appear to be quite consistent 
with the N-body simulations of Bouchet and Hemquist (1992). 

Defining the Fourier transform of the density field, 

6(k) m i f d 3 x6(x)e* x , (17) 

we consider the two- and three-point correlation functions in A:-space, (6(ki)6(k2) ) and ( 6(ki)6{k2)£>(kz ) ), 
which are the Fourier transforms of the spatial two- and three-point functions &( 2 i, 22 ) and £ 3(21 - x 2, 23 ). 
By homogeneity and isotropy, the fc-space moments are non-zero only for 5Z = 0. 

( 6(ki)6(k2) ) = ^ki+kj.O^^l) > ( b{ki)6(k2)S(k$) ) = £k] +k2+kj,0-®(^l * ^"2* ^ 3 ) • (18) 

This defines the power spectrum P(fc) = ( |5(A:) | 2 ) and the bispectrum £123 = B(k\, k2 , £ 3 ). 

Early observations of clustering on small scales (Groth and Peebles 1977) suggested that the galaxy 
two- and three-point functions obey a scaling hierarchy, 

£ 3 ( 0 ?!, X 2 , 23 ) = Q [£ 2 ( 21 , 22 ) 62 ( 252 * * 3 ) + (1 «-> 2) + (2 «-> 3)] (19) 


with Q = constant ~ 1 , roughly independent of the size and shape of the triangle formed by the points 
X\, X21 23 . If the scaling of equation (19) holds exactly, then the hierarchical 3-point amplitude Q is 
also related to the bispectrum by the fc-space version of equation (19) (Fry and Seldner 1982), 


Q = 


B 123 

Pi P2 + Pi P3 + P2P3 


( 20 ) 


with P, = P{ki). We will consider equation (20) as the definition of the amplitude Q, even if it is not 
constant. In the strongly non-linear regime 6 1 , N-body simulations of CDM and power-law spectrum 
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models do seem to display the approximate shape- and size-indepedence of equation (19) (Fry, Melott, 
and Shandarin 1993). However, in second-order perturbation theory in the mildly non-linear regime, 
while Q as defined in equation (20) obeys the scaling with size, it does depend on the shape of the 
configuration in fc-space. 

To calculate the three-point function in the weakly non-linear regime, one expands the pertur- 
bation equations in powers of 6 , 6 (x,t) = 6 { 1 ) {x,t) + M 2 >(x,t) + where is the linear so- 
lution, and <5 (2 > = C?(6 (1) ) 2 is the second-order solution, obtained by using the linear solution in 
the source terms. For Gaussian initial fluctuations, the three-point function vanishes to linear or- 
der, (^ 1 ^(xi)M 1 ^(X2)^ 1 ^(*3) ) = 0, and the lowest order contribution to the bispectrum is £i 23 = 
(^^(fci)^ 1 ^^)^ 2 ^^) ) +(1 *-* 3) + (2 *-> 3), with the result 


Bl23 ~ 


'10 /fcl-fc 2 \ 

(b. + h\ + 4 1 

f ^1 • &2\ 2 

_ 7 + V k x k 2 ) 

U + J + ?i 

l / 


P 1 P 2 + (1~3)+ (2 <-► 3) 


( 21 ) 


(Fry 1984). Strictly speaking, this result holds for initially Gaussian fluctuations in a matter-dominated 
universe with fi = 1, but the work of Juszkiewicz and Bouchet (1991) shows that the dependence of 
the three-point function on fl is extremely slight. A particular case of importance is that of equilateral 
triangle configurations in fc-space, k x = k 2 = *3, for which Q(k) = Q& = 4/7, independent of P(k). 
The independence of this result of the power spectrum makes it a useful quantity for distinguishing 
gravitational from non-gravitational (e.g., bias) effects. (In general, for other configurations or averages 
over configurations, there will be a small dependence on P(k).) In section IV, we will see how this result 
is modified by constant and scale-dependent bias, and compare these predictions with observations. 

Another useful and increasingly popular characterization of the three-point amplitude, which does 
depend on P(k ), is the hierarchical averaged amplitude S 3 , 


Sz(V) = 


MV) = ( 6 3 ( x ; V) ) 
W) ( 6 Hx-,V )) 2 ' 


( 22 ) 


Here f 2 and £3 are the 2-point and 3-point density correlation functions averaged over a window function 
W(r) of characteristic volume V: 


MV) = ^ j J d z r x d 3 r 2 MV\-r 2 \)W{r x )W(r 2 ) 

MV) = J J J d 3 r x d 3 r 2 d 3 r 3 ^(r x ,r 2 ,r 3 )W(r x )W(r 2 )W(r 3 ) (23) 

In comparing with model predictions, it is useful to think of 53 as the ratio of moments of the density 
field 6(x; V ) smoothed over the volume V (Cf. equation (22)), 

6 (x;V) = ^ J d 3 r 6 (x + r)W(r) . (24) 

Thus, is just the variance of the smoothed density field, given by equation (14), and £ 3 (V) is 

its skewness. (The smoothing discussed here should not be confused with the smoothed density field 
introduced in the CGF model of equation (5); in the CGF model, the smoothing radius is associated 
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with the physical scale of threshold modulation effects, while here it merely defines the resolution with 
which one observationally probes the density field.) 

Following standard practice, we evaluate 53 with a top-hat window: for the volume V = 47riZ 3 /3, 
W(r) = 1 for r < R and vanishes for r > R\ its Fourier transform W(kR) is given by equation (15). 
In this case, £ 2 and £3 related to the moments of counts in cells of volume V , and the skewness is 
given by 

(S 3 (R)) = jgp J J <* 3 M 3 *2B(*i,*2,|ki +k 2 \)W{k 1 R)W(k 2 R)W(\k 1 +k 2 |i?) . (25) 

In Fig. 2, we plot 53 as a function of the top-hat smoothing radius R for CDM power spectra with 
Qh = 0.5 and Qh = 0.2 (Cf. (13)), using the second-order perturbation theory result (21) for the 
bispectrum (and assuming that the smoothing radius R is much larger than the galaxy smoothing 
radius R g ~ 1 h~ l Mpc). In computing 53 for the low-density model, we have ignored the tiny correction 
for Q ^ 1 (Juszkiewicz and Bouchet 1991). The result for the CGF model, also shown here, will be 
discussed below in section IV. So far as we are aware, these numerical results for 53 for CDM are new. 
(Goroff et al. 1986 roughly integrated S 3 for CDM with a Gaussian smoothing window using Monte 
Carlo integration, and we have also studied S 3 for Gaussian smoothing. Top hat smoothing requires a 
more accurate numerical integrator, and we have checked our integration code by comparing with the 
analytic results of Juszkiewicz and Bouchet (1991) for S 3 for power law spectra-see below). Where our 
results overlap with the N-body results of Bouchet and Hemquist (1992), the agreement is quite good. 
We see that S 3 does vary with scale R in a manner that depends on the shape of the power spectrum, 
because the CDM spectrum is not exactly scale-free. For a scale-free, power-law spectrum P(k) oc fc n , R 
c an be scaled out of the expression for S 3 , i.e., S 3 is a constant, and its value can be found analytically, 
Ss(R) = 34/7 - (n -I- 3) (Juszkiewicz and Bouchet 1991). On the other hand, for a purely unsmoothed 
field, R = 0, W(kR) = 1, the normalized skewness is S3(0) = 34/7, independent of the power spectrum 
(Peebles 1980). 

The hierarchical behavior of the three-point function in perturbation theory extends to higher order 
correlations, so one cam define higher order hierachical amplitudes Qj ~ £j /£ 2 -1 or = O/f 2 1 
which have characteristic amplitudes set by gravitational instability (see Peebles 1980, Fry 1984, Goroff 
et al. 1986, Bemardeau 1992). 

4 Scale-dependent bias and the 3-point correlations 

We now turn to study how the J-point correlation amplitudes, and in particular the three-point function, 
are affected by constant and scale-dependent biasing. Because we consider the 3 -point function £ 3 , 
we must extend the CGF model to the case in which the matter distribution is not just Gaussian 
but hierarchical, i.e., we consider the contribution of second-order gravitational evolution. Fry and 
Gaztanaga (1993a) have shown that the first-order contribution of biasing to £3 is comparable to the 
contribution from second-order gravitational evolution and, thus, it is not consistent to assume a purely 
Gaussian density field. We first consider how the non-local cooperative modulation of the density field 
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affects the 3-point function, and then study how it is further affected by linear and non-linear bias, that 
is, we consider the sequence of transformations 6 — ► 6' — » 6 g . 

4.1 Cooperative bias 

Consider the effect on the 3-point amplitude of the non-local, cooperative linear transformation of the 
density field given in equation (10). The bispectrum of the cooperative field S'(x) is 

B'\2Z = -®123 (1 + *£i)(l + k{?2)(1 + kQz) , (26) 


where Q x = G(k t ) is given by equation (11). The hierarchical 3-point amplitude Q' of the field S', defined 
in equation (20), can be expressed in terms of the 3-point amplitude Q for the underlying density field, 

6: 

Ql _Q {P\P\ I + P\Pj + PjPti (1 + *ffl)(l + «fe)(l + «fo) 

W V PiP 2 (1 + ^i) 2 (1 + «^2) 2 + (1~3) + (2~3) ‘ (27) 

Note that the ratio Q' /Q has no explicit angular dependence in k- space, i.e., it depends only on the 
magnitudes Aq, k 2 , fc 3 . Using this property, we can point to several important limiting behaviors of 
Q' IQ. For example, on small length scales, Aq, k 2 , k 3 > we obviously retrieve Q' = Q, and in the 
opposite limit of large scales (small triangles in fc-space), ki, k 2 , k 3 « Rj 1 , we have Q'/Q ~ 1/(1 + k), 
independent of the power spectrum and the triangle configuration. The other limiting case of interest is 
a triangle with two large sides and one small side, e.g., Aq, k 2 > .R7 1 , k 3 < Rj l : if the power spectrum 
is approximately a power law, P(k) oc k n , then for n > 0 (and A: 3 /Aq( 2 ) < (1 + k )~ 2/n ), Q' IQ *1 + k; 
for n = 0, Q'/Q = 3(1 -I- k)/\ 1 + 2(1 + /c) 2 ]; and for n < 0, Q'/Q 1/(1 + k). 

As noted in section III, an important class of configurations is equilateral triangles in fc-space, 
fci = &2 = ^3 = k, for which 

* -n (1 + *S) 3 Qa 


Q'a=Q& 


(l + Kfil ) 4 (l + «e) ' 


With the Gaussian CGF smoothing window, Q(k) = for scales larger than R s , kR a -C 1, we 

have Q’ a = Q& (1 + k)~ 1 , whereas for scales smaller than R a , kR, » 1 , we have Q' A = Q& = 4/7. For 
the preferred parameter values considered by BCFW to match the APM data, R a = 20 h -1 Mpc and 
k = 2.29, we see that within the range of the weakly non-linear regime, A: -1 ~ 10h _1 Mpc, there is a 
sharp transition from Q' A — Q a to Q' A = 0.3Q^. We explore the observational consequences of this 
behavior in the next section (see Fig. 3). 

It is also of interest to study the normalized skewness of the smoothed cooperative density field, 
S3 (R) = ( S' 3 (x; R) ) / ( S r3 (x\ R) ) 2 , where the Fourier-transform of the top-hat-smoothed cooperative 
density field is 6'(k;R) = W(A:ii)6(fc)[l + K^(A:iJ s )], with W(kR) given by (15) and Q(kR a ) given by 
(11). The function 5 3 (il) is shown, for the CDM Clh = 0.5 spectrum with the canonical CGF parameters 
k = 2.29, R a = 20 h -1 Mpc, as the curve labelled CGF in Fig. 2. As expected, in this case S3 has a 
steeper dependence on R for scales R < R a than either the standard or low-density CDM models, due 
to the rather sharp, non-power-law feature in S'(k) arising from cooperative effects. These different 
behaviors are compared with data in Fig. 4 below. 
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One can also define higher order amplitues, Qj, by { 6 J (k ) ) ~ QjP J l , with Qz = Q. From the 
above arguments it is straightforward to show that, in general, for regular J-sided polygons in fc-space, 

Qj 


Q'j = 


(1 + kQ) j 


-2 


(29) 


Again, for a Gaussian filter Q, we have Q'j — Qj (1 + k)~ j+ 2 for kR g 1 and Q'j = Qj for kR s 1. 
Thus if the underlying density field 6(x) is hierarchical, with Qj approximately constant as a function 
of k, the new field S' is also hierarchical, ( S‘ J ) = Q'j (S' 2 ) J ~ X , with the hierarchical amplitudes Q'j 
varying with scale from Q'j = Qj to Q'j = Qj (1 + k )~ j+ 2 . 


4.2 Non-linear, local bias 

In the previous subsection, we considered the three-point function for the cooperative density field 
6'(x) defined in equation (5). We now want to relate this to the three-point function of the galaxy field 
6 g (x), defined by the arbitrary, local, non-linear transformation of the cooperative field in equation (6). 
Fry k Gaztanaga (1993a) have shown that, in the weakly non-linear limit { 6 2 ) < 1, the hierarchical 
relation between the moments of the density field, (& ) oc ( 6 2 ) J_1 , is preserved under an arbitrary 
local transformation of this form. Nevertheless, the higher order moments of the galaxy field will differ 
quantitatively from the hierarchical amplitudes of the cooperative field. The analysis of Fry k Gaztanaga 
(1993a) is valid as long as the amplitudes of the original field (here, the cooperative field S'(x)) are of 
zeroth order in the two-point function £2 = (^ 2 ), i.e., under the assumption that Q'j = O { S' 2 )°; in 
particular, their results apply even if the original field is not hierarchical in the strict sense that Q'j is 
constant. 

Let the hierarchical amplitudes of the smoothed galaxy field S g (x) = f(S') be denoted by Q g j To 
consider the 3-point galaxy amplitude, Q g = Q g< 3, we must keep terms up to quadratic order in S' in 
the expansion (7) of the biasing function f(S'). Applying the results of Fry and Gaztanaga (1993a), we 
find 

Q g = b~ l (Q' + C2) + 0(S n ) , (30) 

where C 2 = b^/b and b = b\ in equation (7), and Q * , the 3-point amplitude for the cooperative field 
is related to the 3-point amplitude of the underlying density field by equation (27). For example, for 
the high peaks model, in the limit v 1 and 1, the bias function f{6') is exponential, and we 
have C2 = b. This suggests that the C2 term in (30) is of the same order as the Q' term, i.e., that the 
contribution of non-linear bias to the galaxy 3-point function may be comparable to the second-order 
gravitational contribution. For equilateral triangles in fc-space, we can use (28) and (30) to relate the 
galaxy 3-point amplitude Q g ,& to that of the underlying density field, Q&, 

«»'* = r ‘(l%s +C2 ) ' 1311 

with Q a = 4/7. On small scales, fc/vz, fc >• Rj 1 , where cooperative effects are negligible, but 
still in the mildly non-linear regime ( 6 2 ) < 1 , we have Q g ,ss — b~ l (Q + C2), just the result in the 
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absence of cooperative effects. On large scales, k R a 1 , the galaxy 3-point amplitude is Qg.LS — 
6 -1 [(l + k)~ x Q -I- C2]. The fractional change in Q g between large (kR a < 1) and small {kR a > 1) scales 
is thus 


Q S 


\l + kJ bQg 


(32) 


For the case of purely linear biasing, C2 = 0, this gives AQ g /Q g ,ss — */(l + «), independent of b or Q. 

A similar expression can be derived for the hierarchical amplitudes Sj = Zj/ti °f the volume- 
averaged correlation functions. Following the arguments above, the small- to-large-sc ale variation in the 
galaxy 3-point amplitude S g = S 9 , 3 is related to the density amplitude 5 = S3 by 


AS g ( k \ S 

Sg ~\l+KjbS g • (33) 

Again for purely linear bias, this gives A S g /S g ^ss — /c/(l + ac). 

Finally, for the non-CGF models, note that equation (30) relates the galaxy and matter density 
3-point amplitudes with the replacement Q> — > Q (Fry and Gaztahaga 1993a); we will make use of 
this in comparing the CDM models to observations below. It is also worth reiterating that all of our 
results apply to models with initially Gaussian fluctuations. For non-Gaussian models, there is an 
additional first-order contribution to the 3-point amplitude, which can be thought of as a (possibly 
scale-dependent) contribution to the parameter C2. 


5 Comparison with observations 

We now compare the model predictions for the three-point amplitudes with observations from the 
CfA, SSRS, and Perseus- Pisces redshift surveys. As above, we focus on three models: standard CDM 
{£lh = 0.5), low-density CDM (Qh = 0.2), and CGF-modified standard CDM, and we employ the results 
of second-order perturbation theory. Comparison with the observed galaxy amplitudes, S g ^ and Q 9 ,a, 
can in principle be used to constrain the bias parameters b and C2 as well as the CGF parameters k 
and R a . In combination with other observations, e.g., of the galaxy power spectrum, and of the large- 
angle microwave anisotropy as seen by COBE DMR and other experiments, these results can help point 
toward preferred models for large-scale structure. 

5.1 Limits from Qa 

Baumgart and Fry (1991) have estimated the galaxy power spectrum and the Fourier-space three-point 
amplitude for equilateral triangle configurations, Q&, using data from the Center for Astrophysics and 
Perseus-Pisces redshift surveys. It is worth noting that the power spectrum P(k) for these samples 
does show evidence for the extra large-scale power inferred in other spectroscopic (e.g., IRAS) and 
photometric (e.g., APM) surveys. Their results for Q^(fc), averaged over 3 subsamples each from the 
CfA and Perseus-Pisces surveys, are shown in Fig. 3. The errors bars in each bin indicate the variance 
between subsamples, and we only show results for values of the wavenumber away from the strongly 
non-linear regime. 
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The striking feature of these results is the relative constancy of the three-point amplitude over more 
than a decade in wavenumber, k — 0.1 - 1.6 (h~ l Mpc) -1 . Moreover, the observed amplitude of Qa 
over this range is apparently in reasonable agreement with the prediction of second-order perturbation 
theory without cooperative effects, and under the assumption of no bias, b — 1, C 2 = 0, namely Qa = 4/7 
(shown as the short-dash line in Fig. 3). Turning this around, using the perturbation theory relation 
Q g A = b~ l [( 4/7) + c 2 ], one can in principle use the results in Fig. 3 to constrain the parameter space 
of 6 — c 2 for any model with scale-independent bias. In practice, however, the derived constraints 
are not terribly stringent. First, searching this two-dimensional space and treating the data points 
as independent, one finds a minimum x 2 — 25.5 (for 12 data points and a 2-parameter fit, i.e., 10 
degrees of freedom) for c 2 ~ 0.526 - (4/7). This is consistent with a mean value of Qa ~ 0.52 over 
the plotted range of k , close to the expected perturbation theory result of 4/7 = 0.57. In particular, 
for purely linear bias, c 2 = 0, the best fit value of the bias parameter is 6 = 1.1 ± 0.1, consistent with 
the visual impression from Fig. 3. On the other hand, this constraint on the bias parameter space 
should be interpreted with a great deal of caution, since the best fit curve for perturbation theory has 
a chi-squared of 2.5 per degree of freedom, more than 3-<7 above the expected value. A better fit to the 
data would be obtained with a model in which Qa falls gently with increasing k. However, given the 
likelihood that the true data errors are larger than those shown here, it would certainly be premature 
to exclude the perturbation theory result on this basis. 

The statements above apply for local, non-cooperative bias models. On the other hand, as noted in 
section 4.1, the CGF model predicts a dramatic scale-dependence of Q&(k) around the scale kR 3 ~ 1. 
This behavior is shown in Fig. 3 for the 3 CGF parameter choices considered by BCFW, k,R 9 = 
0.84, 10 /i" 1 Mpc (dot-long dash curve), 2.29, 20/i _1 Mpc (solid curve), and 4.48, 30 h~ l Mpc (dot-short 
dash curve). As above, these models are plotted for c 2 = 0, 6 = 1. The ‘smoking gun’ of these models is 
the sharp downturn in Qa on large scales. Since, within the observational errors, no such downturn is 
observed, one can use this to constrain the CGF parameter space. In particular, for R s = 10h~ l Mpc, 
k = 0.84, the CGF model is always a significantly poorer fit to the data than the scale-independent 
bias models. For this choice of CGF parameters, the requirement of a fit that is within l-o of the scale- 
independent models (i.e., a fit with x 2 < 3 per degree of freedom) necessitates a linear bias parameter 
6 > 2.6 and a significant non-linear bias, c 2 > 0.8. In this case, the large linear bias factor suppresses 
the gravitational and cooperative contribution to Q, and the match with the observations is obtained 
chiefly by the non-linear bias. This would make the apparent agreement between the observed Q a 
and the perturbation theory prediction of 4/7 purely coincidental. This behavior is an instance of our 
general conclusion that models with sharply varying scale-dependent bias are forced to uncomfortably 
large values of the linear bias 6. On the other hand, for larger values of the CGF ‘scale of influence 1 R s , 
the 3-point data do not extend to large enough scales for the downturn to be significant. Consequently, 
the R 8 = 20 h~ l Mpc CGF model, when fitted to the Qa(^) data, occupies the same region in the 
two-dimensional 6 — c 2 bias parameter space, with only a slightly higher x 2 than the non-CGF models. 
Clearly, to more strongly constrain or rule out the CGF model, it would be useful to have data on 
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Qa(^) which extends down to k < 0.05 h Mpc -1 ; this should be feasible with currently available 
redshift samples drawn from the IRAS catalog. 

5.2 Limits on S 3 

To compare model predictions to observations of the volume-averaged normalized skewness S3, we use 
the results of the S3 analysis by Gaztanaga (1992) for samples in the CfA and SSRS redshift catalogs 
(we use the largest samples, denoted SSRS115 and CfA92 in Gatanaga 1992). The average over these 
samples is shown in Fig. 4, where we plot S ff , 3 as a function of top-hat smoothing radius (or cell size) 
R. Each data point in Fig. 4 is am average over bins that correspond to different degrees of freedom: 
for a given value of R, the average number of galaxies in that cell size is at least one unit larger than 
in the cell of the next smallest value of R shown in the figure. The errorbars shown here are the larger 
of the intersample dispersion in the given ii-bin and the intrinsic errors in the original samples. From 
Fig. 4, it is apparent that S 9 i 3 (R) ~ 2 is quite constant over the range of R shown, with a variation of 
about 25%. We also show in Fig. 4 the same model predictions for S3 as in Fig. 2, again for the bias 
parameters 6 = 1 and C2 = 0. The reader should mentally note that the curves in Fig. 4 can be shifted 
vertically, and have their slopes magnified or depressed, by changing the values of 6 and C2. 

In Fig. 5 we plot the contours of x 2 for the comparison between the three models and the S 3 
observations in the 6 - C2 parameter space. The 3 contours correspond to x 2 = 5, 8, and 14 for 11 
data points fit with 2 parameters (9 degrees of freedom). Again, because of the way error bars have 
been assigned to the data, we caution against absolute interpretations of these x 2 values; however, the 
difference in x 2 values for different models should provide a measure of the relative goodness of fit to 
the data. In this sense, both CDM models give comparably good fits to the data for values of the bias 
parameter above 6=1, with 6 > 1.8 for the best fit (lowest x 2 contour) range. For a given value of 
6, the non-linear bias C2 is slightly larger for the Qh = 0.2 case than for standard fl/i = 0.5 CDM. For 
the CGF model, on the other hand, the linear bias parameter must satisfy 6 > 2 for a reasonable fit, 
while the best fit range requires 6 > 3. The large value of bias for the CGF model inferred from S3 
is qualitatively similar to the result above from the Fourier-amplitude Q A : fits to the data with large 
values of 6 are in a sense ad hoc, because the agreement is obtained by depressing the gravitational 
contribution and then fitting with the non-linear bias C2 alone. In particular, for C2/6 ~ 0.6, the galaxy 
amplitude S g< 3 ~ 2 ~ 3Q is completely produced by non-linear bias, not by gravitational or cooperative 
effects. Therefore the fit for the CGF model, for which cy/b — 0.6, does not really reflect agreement 
between the data and the CGF model, but rather the possibility that, in any model, the observed signal 
comes from the non-linear component of biasing. 

At this point, it is worth noting several features of the S3 observations. The skewness has been 
measured from other redshift and angular catalogs in addition to those used above; a useful compendium 
of results in the literature is given in Fry and Gaztanaga (1993b). Except for the Lick catalog, the values 
of S3 inferred from other surveys are broadly consistent with those shown in Fig. 4 (e.g., Bouchet, etal. 
1991, 1993, Meiksin, etal. 1992). A second issue concerns redshift distortions of the higher order 
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moments. It is well known that peculiar velocities distort the galaxy power spectrum (Kaiser 1987), so 
that the power measured in a redshift catalog does not precisely represent the clustering power in real 
space. The transformation from the real space to redshift space power spectrum depends on the ratio 
O 0 6 /6. The extent to which this affects higher moments has been somewhat controversial: in N-body 
simulations, Lahav etal. (1993) find that S 3 is significantly distorted in redshift space in the strongly 
non-linear regime, while Coles, etal. (1993) do not see this affect. In their analysis of higher moments 
in the CfA, SSRS, and IRAS 1.9 Jy catalogs, Fry and Gaztanaga (1993b) find that the volume-average 
3-point function £3 is affected by redshift distortions, but that the normalized skewness 53 is insensitive 
to them. This empirical insensitivity justifies our comparison of the model results to the 53 data in 
redshift space. 

We finish this section with some comments about the implications of the Q and S 3 observations for 
the bias parameter(s) and how these compare with other data on large-scale structure. We will focus 
on the CDM models (as opposed to the CGF model). First, as noted above, the Qa observations do 
not significantly constrain the bias parameter b once one allows for non-linear bias (although they do 
imply a relation between 6 and C 2 ). On the other hand, the S 3 observations do appear to favor larger 
values of the bias, b > 1.8, for both CDM models. In a simple bias prescription, for CfA galaxies we 
would expect bo% ~ 1, so that, taken at face value, this constraint on b would imply a low normalization 
amplitude for the CDM models, 03 < 0.56. For standard Qh = 0.5 CDM, this is uncomfortably low 
compared to the amplitude inferred from the COBE DMR measurement of the large-angle microwave 
anisotropy, a%^mr ~ 1 For the low-density CDM model, the <rg amplitude inferred from COBE has a 
large range, depending on the choice of Q and h (Efstathiou, Bond, and White 1992). For example, 
for the choice Q = 0.3, h = 2/3, Efstathiou, Bond, and White (1992) infer <78 ~ 0.7 from COBE, 
closer to the range implied by the S 3 observations. (On the other hand, for lower Q and larger h, e.g., 
Q = 0.2 and h = 1, the COBE value for <7g becomes larger than unity, which is disfavored by the 
3-point data.) While it is tempting to draw conclusions about the viability of different models from 
this comparison, in particular, to argue against standard COBE-normalized CDM, there are potential 
pitfalls which mitigate against making high confidence- level statements of this type. In particular, if 
one focused only on the 53 data in Fig.4 at large R (where the perturbation result is more trustworthy), 
one would conclude that standard Qh = 0.5 CDM fits the data well with b ~ 1, c 2 ~ 0, in agreement 
with the COBE normalization. A conclusion which can be drawn with more confidence from Fig. 5 is 
that the high peaks model prediction c 2 /6 = 1 is inconsistent with the S 3 data for any of the Gaussian 
models we have studied. 

6 Conclusion 

We have studied the three-point galaxy correlations in models of large-scale structure, focusing on the 
CDM model and its variants with extra large-scale power, working in the context of biased galaxy 
formation in second order perturbation theory. In the non-local bias scheme, galaxies form and light up 
in just such a way as to create the illusion of extra power. We have shown that models with effective 
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scale-dependent (or non-local) bias, such as the CGF model, can display the same enhanced large- 
scale power as other variations of standard CDM, but that they break the scaling hierarchy between 
the two- and three-point functions that arises from gravitational evolution. The resulting step in the 
Fourier-space three-point function Qa(£) at the scale k ~ R~ l of the bend in the bias function (which 
produces the extra large-scale power) should provide a strong observational test of scale-dependent 
bias models. However, this step can be partially masked if b is large and if there is significant non- 
linear bias. Consequently, using data currently available, we have shown that the scale-dependent bias 
explanation of large-scale power requires a larger value of the linear bias factor b than in the standard 
CDM model, said a substantial non-lineax bias, in order to account for the observed flatness of the 
three-point amplitudes. 

On the other hand, the three-point amplitudes S 3 and Q do not strongly discriminate between stan- 
dard and low-density CDM with scale-independent bias; this conclusion also extends to the tilted CDM 
and mixed dark matter models. In these cases, however, the 53 data tentatively point to moderately 
large values of the bias, b > 1.8, but more data on large scales is needed to confirm this. We emphasize 
that it is useful to have observational tests using both 53 and Q&, since the former depends on the 
power spectrum while the latter does not. 

For completeness, we note that the CGF and other non-local bias models have other hurdles to 
overcome in addition to the higher moments. In the CGF and related models, the effective bias fac- 
tor increases with lengthscale. On the other hand, recent N-body simulations of CDM incorporating 
hydrodynamics suggest that the bias factor b(k) decreases with lengthscale (Cf. Katz, Hernquist, and 
Weinberg 1992, Fig. 2 of Cen and Ostriker 1992). In addition, the modifications introduced by CGF 
do not apparently address the difficulties which CDM faces with excessive pairwise velocities on small 
scales (Gelb and Bertschinger 1993 and references therein). On the other hand, it would be interesting 
to study whether there might be a cooperative analogue for velocity bias (Couchman and Carlberg 
1992). 
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Figure Captions 


Fig. 1. The two-point spatial correlation function £(r)/(ba s) 2 in linear theory for standard CDM 
(flh = 0.5), low-density CDM (fl/i = 0.2), and CGF-modified standard CDM with k = 2.29, R s = 
20/i _1 Mpc. 

Fig. 2. The volume- averaged normalized skewness 5a(il) in second-order theory is shown as a function 
of top-hat smoothing radius R for the three models of Fig. 1. 

Fig. 3. The Fourier-space 3-point amplitude for equilateral triangles Q&(1 c) is shown as a function of 
wavenumber k. The data points (from Baumgart and Fry 1991) are an average over subsamples from 
the CfA and Perseus- Pisces surveys. The model points are for standard perturbation theory (short 
dashed line, Q& = 4/7), and for the three CGF models discussed by BCFW: k, R s = 0.84, 10 h -1 Mpc 
(dot-long dash), 2.29,20 /i -1 Mpc (solid), and 4.48, 30 /i — 1 Mpc (dot-short dash). For the models, we 
have taken b = 1, C 2 = 0. 

Fig. 4. The volume-average skewness Sz(R) for the same models as in Fig. 2 are shown in comparison 
with data from the CfA and SSRS surveys (from Gaztahaga 1992). The models are shown with 6=1, 

C2 = 0. 

Fig. 5. Contours of x 2 = 5, 8, and 14 (for 9 degrees of freedom) in the b — C 2 parameter space for fits 
of the 3 models to the data in Fig. 4. The darker regions correspond to lower x 2 - ( a ) CDM Qh = 0.5, 
(b) CDM ft/i = 0.2, (c) CGF k = 2.29, R, = 20h _1 Mpc. 
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